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Abstract 

The production of jets is studied in collisions of virtual photons, 7*p and 
7*7*, specifically for applications at HERA and LEP2. Photon flux factors 
are convoluted with matrix elements involving either direct or resolved pho- 
tons and, for the latter, with parton distributions of the photon. Special 
emphasis is put on the range of uncertainty in the modeling of the resolved 
component. The resulting model is compared with existing data. 



1 To appear in the Proceedings of the International Conference on the Structure and Interactions of 
the Photon; Photon 99, 23-27 May 1999, Freiburg im Breisgau, Germany. 



1 Introduction 



The photon is a complicated object to describe. In the DIS region, i.e. when it is 
very virtual, it can be considered as devoid of any internal structure, at least to first 
approximation. In the other extreme, the total cross section for real photons is dominated 
by the resolved component of the wave function, where the photon has fluctuated into a 
qq state. The nature of this resolved component is still not well understood, especially 
not the way in which it dies out with increasing photon virtuality. This dampening is 
likely not to be a simple function of virtuality, but to depend on the physics observable 
being studied, i.e. on the combination of subprocesses singled out. 

Since our current understanding of QCD does not allow complete predictability, one 
sensible approach is to base ourselves on QCD-motivated models, where a plausible range 
of uncertainty can be explored. Hopefully comparisons with data may then help constrain 
the correct behaviour. The ultimate goal therefore clearly is to have a testable model for all 
aspects of the physics of 7*p and 7*7* collisions. As a stepping stone towards constructing 
such a framework, in this paper we explore the physics associated with the production 
of 'high-px' jets in the collision. That is, we here avoid the processes that only produce 
activity along the 7*p or 7*7* collision axis. For resolved photons this corresponds to the 
'soft' or 'low-p_i_' events of the hadronic physics analogy, for direct ones to the lowest-order 
DIS process 7*q — > q. 

The processes that we will study here instead can be exemplified by 7*7* — > qq 
(direct), 7*g — > qq (single-resolved for 7*7*, direct for 7*p) and gg — > qq (double-resolved 
for 7*7*, (single-)resolved for 7*p), where the gluons come from the parton content of a 
resolved virtual photon or from the proton. Note that these are multi-scale processes, 
at least involving the virtuality Q\ of either photon (i = 1, 2) and the p 2 ± of the hard 
subprocess. For a resolved photon, the relative transverse momentum k± of the initial 
7* — > qq branching provides a further scale, at least in our framework. 

Almost real photons allow long-lived 7* — > qq fluctuations, that then take on the 
properties of non-perturbative hadronic states, specifically of vector mesons such as the 
p°. It is therefore that an effective description in terms of parton distributions becomes 
necessary. Hence the resolved component of the photon, as opposed to the direct one. 
That such a subdivision is more than a technical construct is excellently illustrated by 
the x° hs plots from HERA Q. 

The resolved photon can be further subdivided into low-virtuality fluctuations, which 
then are of a nonperturbative character and can be represented by a set of vector mesons, 
and high-virtuality ones that are describable by perturbative 7* — > qq branchings. The 
former is called the VMD (vector meson dominance) component and the latter the anoma- 
lous one. The parton distributions of the VMD component are unknown from first prin- 
ciples, and thus have to be based on reasonable ansatze, while the anomalous ones are 
perturbatively predictable. 

The traditional tool for handling such complex issues is the Monte Carlo approach. 
Our starting point is the model for real photons || and the parton distribution param- 
eterizations of real and virtual photons already present in the Pythia j|] generator. 
Several further additions and modifications have been made to model virtual photons, as 
will be described in the following |5|]. 
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2 The Model 



The cross sections for the processes ep — > eX and ee — > eeX can be written as the 
convolutions [B, [7], 0] 

da(ep - eX) /"/" dQ 2 /* /e (y, Q 2 ) d^p - X) (1) 

£=T,L ^ ^ 

and 

da(ee eeX ) =J2 [[[ [dy^Qjdy, dQ 2 /&( yi , Cft)f*} 9 fa, QftMl&ib - X). (2) 

The flux of photons /(?/, Q 2 ) (see y definition below) from the lepton is factorized from the 
subprocess cross sections involving the virtual photon, 7*p — ► X and 7*7* — > X. The sum 
is over the transverse and longitudinal photon polarizations. For ep events, this factorized 
ansatz is perfectly general, so long as azimuthal distributions in the final state are not 
studied in detail. In e + e~ events, it is not a good approximation when the virtualities 
Qi and Q\ of both photons become of the order of the squared invariant mass W 2 of the 
colliding photons 0. 

When Q 2 /W 2 is small, one can derive |I| §, § 



f T / n 2N _ «em / (1 + (1 — y) 2 1 2m 2 y » 



fi/i(y,Q 2 ) - (4) 

The y variable is defined as the lightcone fraction the photon takes of the incoming 
lepton momentum. The lepton scattering angle Qi is related to Q 2 , where the kinematical 
limits on Q 2 are, unless experimental conditions reduce the 6 range, Q 2 min ~ iz^ m e an d 

Qf.max « (1 - V)S. 

Within the allowed region, the phase space is Monte Carlo sampled according to 
{dQ 2 /Q 2 ) (dy/y) dip, with the remaining flux factor combined with the cross section fac- 
tors to give the event weight used for eventual acceptance or rejection. 

The hard-scattering processes are classified according to whether one or both photons 
are resolved. For the direct process 7^7^ — * ff, f some fermion, the cross sections for 
transverse and longitudinal photons are used [11, ||. Remember that the cross section for 



a longitudinal photon vanishes as Q 2 in the limit Q 2 —>■ 0. 

For a resolved photon, the photon virtuality scale is included in the arguments of the 
parton distribution but, in the spirit of the parton model, the virtuality of the parton 
inside the photon is not included in the matrix elements. Neither is the possibility of 
the partons being in longitudinally polarized photons (see below, however). The same 
subprocess cross sections can therefore be used for direct 7*p processes and for single- 
resolved 7*7* ones. Both transverse and longitudinal photons are considered for the QCD 
Compton 7*q — > gq and boson-gluon fusion 7*g — > qq processes |L2], ||. The matrix 
elements are then convoluted with parton distributions. 

Finally we come to resolved processes in 7*p and doubly-resolved ones in 7*7*. There 
are six basic QCD cross sections, qq' — ► qq', qq — ► q'q', qq — > gg, qg — > qg, gg — > gg and 
gg — > qq. The same subprocess cross sections as those known from pp physics [[H| can 
therefore be used. Again, a convolution with parton distributions is necessary. 
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One major element of model dependence enters via the choice of parton distributions 
for a resolved virtual photon. These distributions contain a hadronic component that is 
not perturbatively calculable. It is therefore necessary to parameterize the solution with 
input from experimental data, which mainly is available for (almost) real photons. In 
the following we will use the SaS distributions [Q, which are the ones best suited for our 
formalism. Another set of distributions is provided by GRS ||14|| , while a simpler recipe 



for suppression factors relative to real photons has been proposed by DG |[Tj| 
The SaS distributions for a real photon can be written as 

mx,^) = ^^r a v ^^Ql) + ^J2 2 < f^/rW;* 2 )- (5) 

V JV q J Qo 

Here the sum is over a set of vector mesons V = p°, u, </>, J/i/) according to a vector- 
meson-dominance ansatz for low-virtuality fluctuations of the photon, with experimentally 
determined couplings 4na ern / f v . The higher-virtuality, perturbative, fluctuations are 
represented by an integral over the virtuality k 2 and a sum over quark species. We will 
refer to the first part as the VMD one and the second as the anomalous one. 

From the above ansatz, the extension to a virtual photon is given by the introduction 
of a dipole dampening factor for each component, 

+ r ^ZKj^(^p^,^). (6) 

Thus, with increasing Q 2 , the VMD components die away faster than the anomalous ones, 
and within the latter the low-k 2 ones faster than the high-A; 2 ones. 

Since the probed real photon is purely transverse, the above ansatz does not address 
the issue of parton distributions of the longitudinal virtual photons. One could imagine 
an ansatz based on longitudinally polarized vector mesons, and branchings 7£ — > qq, but 
currently no parameterization exists along these lines. We will therefore content ourselves 
by exploring alternatives based on applying a simple multiplicative factor R to the results 
obtained for a resolved transverse photon. As usual, processes involving longitudinal 
photons should vanish in the limit Q 2 —* 0. To study two extremes, the region with a 
linear rise in Q 2 is defined either by Q 2 < fi 2 or by Q 2 < m 2 , where the former represents 
the perturbative and the latter some non-perturbative scale. Also the high-Q 2 limit is 
not well constrained; we will compare two different alternatives, one with an asymptotic 
fall-off like 1/Q 2 and another which approaches a constant ratio, both with respect to the 
transverse resolved photon. (Since we put f2*(x, fi 2 , Q 2 ) = for Q 2 > fi 2 , the R value 
will actually not be used for large Q 2 , so the choice is not so crucial.) We therefore study 
the alternative ansatze 

R 1 ( y ,Q 2 ,ix 2 ) = i+a , ^ ^ 9 ? T /t r , ; 9 : , m 



4/. 2 Q 2 




(/i 2 + Q 2 ] 


2 / 7 T /^,Q 2 ) 


4Q 2 


fyi(y,Q 2 ) 


(/i 2 + Q 2 ] 


f^Q 2 ) 


4Q 2 


fii(y,Q 2 ) 



R 2 (y,Q 2 ,V 2 ) = 1 

R3iM) = 1+a __ ? _g^ (9) 
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with a = 1 as main contrast to the default a = 0. The y dependence compensates for the 
difference in photon flux between transverse and longitudinal photons. 

Another ambiguity is the choice of fi 2 scale in parton distributions. Based on various 
considerations, we compare three different alternatives: 

2 2 s Ql Q2 2 2 1 r\l 1 /02 2 o 2 / 1n \ 

ni = P± 1 , /X 2 = P± + Ql + <5 2 , ^3 = 2 ^1- ( 10 ) 

Only the second alternative ensures f2*(%, A* 2 , Q 2 ) > for arbitrarily large Q 2 ; in all other 
alternatives the resolved contribution (at fixed p±) vanish above some Q 2 scale. The last 
alternative exploits the well-known freedom of including some multiplicative factor in any 
(leading-order) scale choice. When nothing is mentioned explicitly below, the choice fi 2 
is used. 

The issues discussed above are the main ones that distinguish the description of pro- 
cesses involving virtual photons from those induced by real photons or by hadrons in 
general. In common is the need to consider the buildup of more complicated partonic 
configurations from the lowest-order 'skeletons' defined above, (i) by parton showers, (ii) 
by multiple parton-parton interactions and beam remnants, where applicable, and (Hi) by 
the subsequent transformation of these partons into the observable hadrons. The latter, 
hadronization stage can be described by the standard string fragmentation framework 
HJ, followed by the decays of unstable primary hadrons, and is not further discussed 



here. The parton shower, multiple-interaction and beam-remnant aspects are discussed 
elsewhere H. 



3 Comparisons with Data 

In this section the model is compared with data. We will not make a detailed analysis 
of experimental results but use it to point out model dependences and to constrain some 
model parameters. 



3.1 Inclusive Jet Cross Sections 

Inclusive ep jet cross sections have been measured by the HI collaboration [17| in the 
kinematical range < Q 2 < 49 GeV 2 and 0.3 < y < 0.6. For da ep /dE]_ data is available 
in nine different Q 2 bins; two of them are shown here, Fig. [I], with similar results for 
the intermediate bins. The plots were produced with the HzTool [I8[] package. The E* ± 



and rj* are calculated in the 7*p centre of mass frame where the incident proton direction 
corresponds to positive 77*. The SaS ID parton distribution together with a few different 
fii scales are used to model the resolved photon component. 

In the highest Q 2 bin the direct component gives the dominant contribution. However, 
the resolved component is not negligible. All the scales a*i depend on the photon virtuality. 
This gives a larger resolved component in this region as compared to the conventional 
choice, (jl = p±. In the low Q 2 bin the a*i and a*2 scales do not differ much from p± and 
the cross sections are in nice agreement with data. The cross section with the a*3 scale 
overshoots the data in this region. 

Changing the photon parton distribution from SaS ID to SaS 2D will give a slightly 
lower result for the low-Q 2 bins. A comparison has been made with the GRS LO [f[4j| 
parton distribution U and, in its region of validity, differences are small. One could imag- 
ine larger differences for 7* parton distributions that from the onset are more different. 
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Figure 1: The differential jet cross section da ep /dE]_ for jets with —2.5 < rj* < —0.5 
and 0.3 < y < 0.6. a = indicates that only transversely polarized resolved photons are 
considered. 

Using a parton distribution for a real photon cannot describe the Q 2 dampening in the 
distributions shown in this section. 

Using CTEQ 3L instead of GRV 94 LO (which is default in Pythia) as the proton 
parton distribution reduces the result in some E* ± bins by half. The GRV 94 HO parton 
distribution give a slightly lower result (as compared to GRV 94 LO). The differences 
mainly come from the gluon distributions, that are not yet so well constrained from data. 
In the modeling of the parton distributions, it is a deceptive accident that the more 
well-known proton parton distribution gives a larger uncertainty than the photon one. It 
offers a simple example that also phenomenology of other areas may directly influence the 
interpretation of photon data. 

The OPAL collaboration has measured inclusive one-jet and two-jet cross sections in 
the range \r]i et \ < 1 and requiring E 3 ^ to be larger than 3 GeV | 19|| . The centre of mass 



energies were 130 and 136 GeV. The inclusive jet cross sections as a function of E 3 ± or 
rf et are compared with data 0, with events generated at ^/s^ = 133 GeV. 

At low E*l the double-resolved events are dominating and at larger E 1 ^ it is the direct 
processes since more energy goes into the hard scattering in the latter case. For single- 
resolved events, the SaS ID VMD component dies out much quicker with increasing E 3 ^ 
than the SaS 2D one which is comparable with the direct-anomalous events at high E± . 
For both cases, at high E 3 ^, the direct -anomalous components give the same order of 
magnitude contribution to the cross section as the double-resolved events. The biggest 
difference between the two parton distributions can be seen at low E 3 ^ and for the |^ et | 
distributions, where the double-resolved events dominate; it is a reflection of the difference 
in normalization among the contributions. For the SaS 2D case, this kinematical region 
makes the VMD component more important than the anomalous one; as a consequence 
multiple interactions play an important role. The double-resolved contribution for SaS 2D 
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without multiple interaction is reduced by half. Clearly, for the SaS ID case the opposite 
is true: the importance of the components are reversed. In the region of high E 3 ^ , where 
the direct events dominate, the model is undershooting data. On the other hand, there 
is nice agreement with data for the \7f et \ distribution when using SaS 2D. 



3.2 Forward Jet Cross Sections 

Jet cross sections as a function of Bjorken-x, xbj, for forward jet production (in the proton 
direction) have been measured at HERA pOj . The objective is to probe the dynamics 
of the QCD cascade at small XBj- The forward jet is restricted in polar angle w.r.t. the 
proton and the transverse momenta j?± should be of the same order as the virtuality of 
the photon, suppressing an evolution in transverse momenta. If the jet has a large energy 
fraction of the proton there will be a big difference in x between the jet and the photon 
vertex; XBj <C Xj et , allowing an evolution in x. The above restrictions will not eliminate 
the possibility of having a resolved photon, although the large Q 2 values are not in favour 
of it. 

The HzTool routines were used to obtain the results in Fig. [|. A larger forward jet cross 
section is obtained with a stronger Q 2 dependence for the hard scale, with n\ = p 2 ± + Q 2 in 



best agreement with data ||21|| . The choice of scale does not only affect the resolved photon 
contribution but also the direct photon, arising from the scale dependence in the proton 
parton distribution. The rather large Q 2 values, Q 2 ~ (p 3 ±) 2 , suppresses VMD photons 
and favours the SaS ID distribution which is the one used here, though the difference is 
small. 
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Figure 2: Forward jet cross section as a function of x. The results with three different 
alternatives of longitudinal resolved photons Ri are compared with purely transverse ones, 
a = 0, and data from HI. 

Note that the /z 3 scale undershoots the forward jet cross section data and overshoots 
the inclusive jet distributions at low Q 2 , so it is not a real alternative. As a further 



6 



check, with more data accumulated and analysed, the (p^) 2 /Q 2 interval could be split 
into several subranges which hopefully would help to discriminate between scale choices. 

With the experience of forward jets at HERA, we suggest a similar study at LEP. 
The optimal kinematical and forward jet constraints have to be set by each collaboration 
itself; we will only estimate the order of magnitude for the cross section and point out 
uncertainties in the model. 

Comparing with forward jets at HERA, one of the leptons will play the role of the 
proton. Some of the constraints can be taken over directly, for example, x- ]e t = E- ]et / E e > 
0.035 and 0.5 < {^) 2 /Q 2 < 2. To fulfill the jet selection one of the leptons has to be 
tagged in order to know the virtuality of the photon. To obtain a reasonable number 
of events the other lepton is not tagged. With a centre of mass energy of 200 GeV, the 
smallest accessible xbj = is around 10~ 4 , where Q 2 and y is calculated from the tagged 
electron, omitting the virtuality of the other photon. In a more sophisticated treatment 
also double-tagged events are analyzed; then one of the photons plays the role of a proton 
and the forward jet should be defined with respect to one of the photons. 

As for the case at HERA, the /i2 scale gives the largest forward jet cross section, 
about twice as large as with the /ii scale |5]. Most of the differences arise from the 
double-resolved events. Double-resolved and single-resolved events, where the resolved 
photon give rise to the forward jet, dominate the forward jet cross section. At low x, 
for the /i2 scale, the double-resolved contribution is close to an order of magnitude larger 
than the direct one. For the fii scale it is about a factor of four. As for the case at HERA, 
the rise of the forward jet cross section at small x is dominated by resolved photons. A 
study like this at LEP could be an important cross check for the understanding of resolved 
photons and that of small-x dynamics. 

3.3 Importance of longitudinal resolved photons 

In this section we will study the importance of longitudinal resolved photons. A sensi- 
ble Q 2 ~clependent scale choice, /xi, together with the SaS ID distribution will be used 
throughout. 

With a = 1 the different alternatives are shown in Fig. [3] for the da ep / &E* ± distributions 
together with the result from pure transverse photons, i.e. a = 0. The importance of the 
resolved contributions decreases with increasing Q 2 which makes the asymptotic behaviour 
less crucial. The onset of longitudinal photons governed by the R\ and R2 alternatives 
are favoured whereas the R3 one overshoots data in the context of the other model choices 
made here. 

In Fig. |2| the same alternatives are shown for the forward jet cross sections. With 
this scale choice, fi\, none of the longitudinal resolved components (together with the 
direct contribution) are sufficient to describe the forward jet cross section. The resolved 
contribution with R 3 is about the same as the one obtained with the scale fi 2 = "p\ + Q 2 
(without longitudinal contribution); the difference in the total results originates from the 
difference in the direct contributions. With R\ and a = 1, the fi2 scale (not shown) 
overshoots the data. The above study indicates, as expected, that longitudinal resolved 
photons are important for detailed descriptions of various distributions. It cannot by itself 
explain the forward jet cross section, but may give a significant contribution. Combined 
with other effects, for example, different scale choice, parton distributions and underlying 
events, it could give a reasonable description. The model(s) so far does not take into 
account the difference in x distribution or the k 2 scale (of the 7* — > qq fluctuations) be- 
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tween transverse and longitudinal photons. As long as the distributions under study allow 
a large interval in x the average description may be reasonable. In a more sophisticated 
treatment these aspects have to be considered in more detail. 



4 Summary and Outlook 

The plan here is to have a complete description of the main physics aspects in 7p and 
77 collisions, which will allow important cross checks to test universality of certain model 
assumptions. As a step forward, we have here concentrated on those that are of importance 
for the production of jets by virtual photons, and are absent in the real-photon case. While 
we believe in the basic machinery developed and presented here, we have to acknowledge 
the many unknowns — scale choices, parton distribution sets (also those of the proton), 
longitudinal contributions, underlying events, etc. — that all give non- negligible effects. 
To make a simultaneous detailed tuning of all these aspects was not the aim here, but 
rather to point out model dependences that arise from a virtual photon. 

When Q 2 is not small, naively only the direct component needs to be treated, but in 
practice a rather large contribution arises from resolved photons. For example, for high- 
er 2 studies like forward jet cross sections, Fig. |2], or inclusive differential jet cross sections, 
Fig. [3]. Resolved longitudinal photons are poorly understood and the model presented 
here can be used to estimate their importance and get a reasonable global description. 
Longitudinal effects are in most cases small but of importance for fine-tuning. 

The inclusive 7*7* one-jet and two-jet cross sections are well described except for the 
high Ej* region of the E 3 ^ distribution || . In this region, the direct events are dominating. 
Currently, owing to the lesser flexibility in the modeling of the direct component, we do 
not see any simple way to improve the model. The factorized ansatz made for the photon 
flux is expected to be valid in this kinematical range; interference terms are suppressed 
by QlQl/W 2 *^. However, differences in the application of the cone jet algorithm may 
affect the results. 

The forward jet cross section presented by HI is well described by an ordinary 
parton shower prescription including the possibility of having resolved photons. The 
criteria that the p^* should be of the same order as Q 2 , makes the scale choice crucial and 
A*2 = v\ + Q 2 is favoured by data, as concluded in [|2T| . With this experience we predict 
the forward jet cross section to be obtained at LEP ||. With more data accumulated and 
analyzed, the {j?±) 2 /Q 2 interval could be split into several subranges, which hopefully 
would help to discriminate between different scale choices. 

Multiple interactions for the anomalous component are not yet included, and is not 
expected to be of same importance as in the VMD case. However, for low k 2 fluctuations 
it may be important, especially for SaS ID, and need to be investigated. 

After this study of jet production by virtual photons it is natural to connect it together 
with low-p_i_ events. Clearly, a smooth transition from perturbative to non-perturbative 
physics is required. Further studies are needed and will be presented in a future publica- 
tion. 
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